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SECTION  I 


INTRODUCTION 

This  literature  summary  reviews  both  experimental  and  analytical  studies 
relevant  to  the  Air  Force  Weapons  Laboratory  (AFWL)  Theoretical  Investigation 
of  Failure  Criteria  for  Reinforced  Concrete  Structures.  First  for  the 
experimental,  and  then  for  the  analytical  studies,  the  summary  generally 
begins  with  earlier,  simpler  studies  and  proceeds  to  the  recent,  more  refined 
studies.  A list  of  the  references  used  in  this  literature  survey  concludes 
this  summary. 
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SECTION  II 
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EXPERIMENTAL  STUDIES 

Static  uniaxial  and  biaxial  tests  comprise  most  of  the  experiments  on 
concrete.  Relatively  few  tests  consider  triaxial,  dynamic,  or  reinforcement 
bond  slip  effects. 

1.  STATICALLY  LOADED  BIAXIALLY  STRESSED  PLAIN  CONCRETE 

Many  experimental  studies  have  been  conducted  over  the  last  70 
years  to  confirm  or  propose  new  theories  governing  the  strength  of 
statically  loaded  biaxially  stressed  plain  concrete. 

In  1911,  Karman  (ref.  1)  and  Baker  (ref.  2)  tested  the  validity  of 
Mohr's  failure  theory  for  brittle  materials.  They  concluded  that  (1) 
depending  on  the  stress  combination,  either  a tensile  splitting  or  a 
shearing  and  sliding  action  governs  failure,  and  (2)  contrary  to 
Mohr's  theory,  the  intermediate  principal  stress  magnitude  influenced 
the  failure. 

Brandtzaeg  and  Brown  (ref.  3)  performed  similar  but  more  refined 
tests  in  1929.  They  concluded  that 

a.  The  biaxial  is  at  least  as  great  as  the  uniaxial  compressive 
strength  of  concrete. 

b.  Three  load  phases  occur  in  uniaxial  compression 

(1 ) Linear  elastic, 

(2)  Inelastic  with  increased  deformation  rate  and  increased 
Poisson's  ratio. 


’ 
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(3)  Internal  cracking  at  75%  to  85%  of  the  maximum  load, 
followed  by  a great  increase  in  lateral  deformation  and 
volume  under  continued  loading. 

c.  In  multiaxial  compressive  stress  - ates  frictional  resistance 
against  sliding  motion  acts  on  a greater  surface  area  than  in 


j 
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uniaxial  compression,  which  could  explain  the  increase  in  applied 
stress  required  to  start  splitting  in  the  multiaxial  case, 
d.  Sliding  may  be  a bond  failure  between  the  cement  and  aggregate; 
the  bond  slip  may  be  partly  adhesive  and  partly  frictional. 


These  conclusions  agreed  somewhat  with  the  internal  friction  theory 
of  sliding  failure.  However,  the  breakdown  process  and  resulting  large 
deformations  could  not  be  reconciled  with  Mohr's  theory  or  any  theory  based 
on  the  assumption  of  sliding  on  planes  of  least  resistance. 

Bresler  and  Pister  (refs.  4,  5,  and  6)  studied  concrete  shear-compres- 
sion stress  states  to  verify  their  octahedral  shearing  stress  theory, 
which  assumed  that  the  octahedral  shearing  stress  at  failure  is  a function 
of  the  octahedral  normal  stress  x0  = f(a0).  For  their  range  of  data  the 
relation  was  linear. 

Freudenthal  (ref.  7)  used  a closely  related  theory  based  on  the  effec- 
tive stress  expressed  in  terms  of  the  octahedral  stress. 

McHenry  and  Kami  (ref.  8)  performed  tensile  compressive  biaxial 
tests  and  also  found  a linear  relation  between  octahedral  shearing  and 
octahedral  normal  stress.  Mohr's  theory  of  failure  was  found  not  appli- 
cable. 
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Later  studies  by  Bellamy  (ref.  9)  and  Goode  and  Helmy  (ref.  10) 

•] 

found  failure  criteria  of  the  form  x = f(a  ) to  be  no  more  accurate  than 

o o 

Mohr's  theory, 

Foppl  (ref.  11),  Wastlund  (ref.  12),  Glomb  (ref.  13)  and  Weigler 
and  Becher  (ref.  14)  performed  biaxial  tests  on  concrete  cubes  to  study 
the  strength  and  failure  aspects  of  the  material.  Again,  compressive 
strength  was  found  to  be  greater  in  biaxial  stress  states  than  in  uniaxial 
compression. 

Vile  (ref.  15)  performed  biaxial  compressive  tests  on  concrete  plates 
in  which  the  stress  strain  curves  were  found  to  be  similar  to  the  uniaxial 
stress  strain  curve. 

Robinson  (ref.  16)  continued  Vile's  tests.  From  the  test  results  it 
was  observed  that:  (1)  the  ultimate  strength  of  concrete  is  increased  by 
lateral  loads,  (2)  increase  In  the  lateral  compressive  load  increases  the 
load  at  which  major  microcracks  occur;  however,  this  increase  was  not  as 
extensive  as  that  of  the  Increase  in  the  ultimate  strength  of  concrete, 

(3)  the  loads  at  which  major  microcracks  occur  decrease  with  increasing 
water/cement,  and  the  aggregate/cement  ratios,  in  contrast  to  the  strength 
of  concrete,  providing  there  is  a lateral  compressive  load.  The  effect  was 
not  so  noticeable  for  the  cracking  loads. 

The  work  suimarlzed  above  is  one  of  the  early  investigations  on  the 
microcracking  of  concrete  under  combined  stresses.  However,  the  results 
obtained  from  this  work  and  microcrack  observations  do  not  appear  to  be  ( 

conclusive. 
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was  obtained  at  a principal  stress  ratio  of  o2/aj  = 0.5.  For  equal 
compressive  stresses  o1  = a2  the  strength  Increase  was  about  16 
percent. 

2.  Specimens  loaded  in  uniaxial  compression  showed  numerous  micro- 
cracks parallel  to  the  direction  of  the  applied  load. 

3.  Specimens  subjected  to  biaxial  compression  showed  microcracks 
parallel  to  the  free  surfaces  of  the  specimens. 

4.  External  major  stress  vs.  strain  relations  were  plotted  for 
various  a2  transverse  confining  loads.  Generally,  for  a given 
Oj  compressive  stress,  the  strains  in  the  first  direction  are 
decreasing  with  increasing  values  of  a2  compressive  stress. 

However,  for  equal  compression  Oj  = a2,  and  at  about  ultimate 
load,  larger  strains,  as  compared  to  those  obtained  by  stress 
ratio  a2/a!  = 0.5,  were  recorded.  Figure  1 presents  a plot  of 
the  general  shape  of  the  elastic  and  ultimate  stress  levels  as 
as  defined  by  Kupfer.  Kupfer  more  recently  (ref.  18)  conducted 
further  research  in  an  attempt  to  describe  the  deformations  of 
plain  concrete  under  biaxial  stresses  by  means  of  mathematical 
expression  in  order  to  predict  the  strains  resulting  from  arbi- 
trary plane  stress  states.  Such  Information  is  required  for 
the  formation  of  material  stiffness  matrices  In  finite  element 
computations.  The  effort  makes  use  of  the  decompositions  of 

the  stress  and  strain  states  into  volumetric,  or  hydrostatic  portion, 
and  a distortional or  deviatory  portion.  Figure  2 shows  a 
plot  of  the  octahedral  shear  stress-strain  curves  for  different 
stress  ratios.  Kupfer  concludes  that  since  the  curves  for  dif- 
ferent stress  ratios  coincide  closely,  a unique  relationship  which 


II 


Is  independent  of  the  principal  stress  ratio  exists.  This  assump- 
tion is  commonly  applied  to  the  theory  of  Inelastic  isotropic 
solids.  Kupfer  continues  to  present  plots  and  relationships  relating 
the  octahedral  shear  stress  and  strain  to  the  material  tangent  bulk 
and  shear  moduli.  Figures  3 and  4 present  plots  of  octahedral  shear 
stress  to  shear  and  bulk  moduli.  The  associated  mathematical  relation- 
ships for  the  tangent  terms  are: 
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where  a,  m,  c,  p,  f are  material  constants  and 
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A suntna ry  is  given  In  the  report  for  the  application  of  this 
procedure  to  the  finite  element  approach. 


Buyukozturk  (ref.  19)  made  some  refinement  of  previous  work  for  the 
behavior  of  concrete  under  uniaxial  loading,  but  his  main  emphasis  was  on 
the  short-term  compressive  biaxial  cases.  Within  the  compressive  region 
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of  biaxial  stress  combinations,  four  different  principal  stress  ratios, 

a /o  , were  chosen,  i.e.,  0.0,  0.2,  0.5,  and  1.0.  For  the  experimental 
2 1 

investigation,  a biaxial  load  testing  device  was  built  which  permitted 
independent  application  of  stress  in  two  perpendicular  directions  and 
which  permitted  x-ray  examination  of  the  specimen  for  cracking  during 
test,  without  the  necessity  of  disturbing  the  specimen  or  removing  the 
load.  In  loading,  special  care  was  taken  to  minimize  the  restraint 
offered  by  loading  plates  due  to  deformation  in  the  direction  normal 
to  the  load.  For  this  purpose,  a brush-like  arrangement  originally  pro- 
posed by  Hilsdorf  (ref.  17)  was  used.  In  this  investigation  the  concrete 
was  idealized  as  square,  flat  specimens  consisting  of  nine  circular  discs 
of  aggregate  embedded  in  a mortar  matrix.  Buyukozturk's  specimens  were 
tested,  and  total  deformations  in  two  principal  directions  were  measured. 
Deformation  versus  external  load  relationships  were  established.  The  effect 
of  the  second  principal  stress  on  the  strength  and  deformations  in  the 
first  principal  direction  was  emphasized.  Using  X-ray  procedure,  micro- 
cracks were  examined  during  actual  loading  process.  Local  strains  were 
measured  at  certain  critical  locations  in  the  specimen. 

On  the  basis  of  Buyukozturk's  studies,  certain  important  conclusions 
can  be  drawn  as  follows: 

1.  Significantly  higher  strength  is  attainable  for  a given  material 
in  biaxial  loading  than  in  uniaxial  loading.  The  strength  in- 
crease is  dependent  upon  the  ratio  of  the  principal  stresses, 
and  appears  to  be  a maximum  at  a stress  ratio  of  about  0.5, 
diminishing  somewhat  as  the  ratio  is  increased  to  unity. 


2.  The  stiffness  In  the  first  principal  direction  is  significantly 
increased  by  the  introduction  of  principal  stress  in  the  perpen 
dicular  direction;  i.e.,  compressive  deformation  in  the  first 
direction  is  substantially  reduced  by  compressive  stress  in  the 
second  direction. 

3.  When  spacing  of  aggregate  of  one  size  is  decreased,  uniaxial 
strength  is  decreased  but  biaxial  strength  is  not  affected 
appreciably. 

4.  In  the  uniaxial  case,  failure  occurs  by  progressive  micro- 
cracking, starting  at  the  aggregate-mortar  interfaces,  and 
later  extending  as  tensile  cracks  through  the  mortar.  In 
the  biaxial  case  cracks  perpendicular  to  the  plane  of  the 
specimen  were  neither  observed  by  experiment  nor  predicted 
by  analysis. 

5.  In  the  uniaxial  case,  ultimate  failure  occurs  by  splitting  in 
the  planes  perpendicular  to  the  face  of  the  specimen  and 
parallel  to  the  load.  In  the  biaxial  case,  ultimate  failure 
occurs  by  splitting  along  a plane  parallel  to  the  face  of  the 
specimen. 

Liu  (ref.  20)  extended  Buyukozturk's  work  by  developing  a more  com- 
plete concrete  model  specimen.  Instead  of  using  nine  circular  disks  of 
aggregate,  Liu's  model  contained  thirty  with  three  different  disc  sizes 
Eight  such  models  were  tested  under  short  term  static  compressive  loads, 
applied  from  zero  through  the  full  range  of  biaxial  compressive  loading 
to  ultimate  strength.  The  failure  envelopes  obtained  through  the  tests 
showed  good  agreement  with  those  of  Buyukozturk's.  Liu  also  subjected 
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real  concrete  specimens  to  biaxial  compression.  Eighty- four  specimens 

with  five  different  mixes  were  tested  at  four  stress  ratios,  i.e., 

a /a  = 0.0,  0.2,  0.5,  and  1.0.  Three  of  the  concrete  mixes  were  used 
2 1 

to  study  the  influence  of  water-cement  ratio  and  aggregate-cement  ratio. 

To  study  the  effect  of  coarse  aggregate  (concrete  vs.  mortar),  the  fourth 
specimen  type  was  a mortar  mix  with  the  same  water-cement  ratio  as  one  of 
the  first  three  concrete  mixes.  To  study  the  effect  of  size  of  coarse 
aggregates,  the  fifth  concrete  mix  had  1 inch  maximum  size  of  aggregate. 

The  results  of  testing  these  84  speciments  showed  good  agreement  with 
the  idealized  model  of  concrete  used  previously.  Studying  the  different 
mix  results  showed  that  ultimate  strength  of  mix  five,  the  mortar  mix,  was 
found  to  be  only  8 to  13  percent  lower  than  that  of  concrete.  The  effect 
of  using  larger  aggregate  did  decrease  the  ultimate  strengths,  but  not 
appreciably.  Figure  5 presents  an  ultimate  stress  plot  of  the  loading 
results  of  Liu's  tests  with  the  different  types  of  mixes. 

A review  of  many  stress-strain  relationships  proposed  by  previous 
researchers  in  concrete  under  biaxial  loading  was  conducted  in  this  paper. 

A new  relationship  was  proposed  by  Liu  for  concrete  under  biaxial  compression, 
biaxial  tension,  or  biaxial  tension-compression.  Comparisons  of  values  from 
the  proposed  rule  with  experimental  data  provided  good  agreement.  This 
relation  is  discussed  further  in  the  analytical  portion  of  this  report. 

Until  recently,  little  biaxial  testing  has  been  carried  out  in  which 
tension  was  one  of  the  loadings  because  of  the  difficulty  involved  in 
applying  a tensile  load  to  conrete.  As  a result,  some  of  the  early. attempts 
at  studying  compression-tension  loading  involved  dubious  testing  techniques. 
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Carlno  (ref.  21),  in  his  paper  published  in  1974,  gives  a brief  history 
of  the  testing  problems  and  the  evolutionary  process  to  arrive  at  today's 
state  of  the  art. 

Carino  states  the  following  conclusions  from  his  review  of  previous 
direct  tension  tests  on  concrete: 

1.  Concrete  possesses  a discontinuity  point  in  direct  tension, 
which  occurs  at  about  70  percent  of  the  short-term  ultimate 
strength. 

2.  The  point  of  discontinuity  corresponds  to  the  formation  of  the 
first  crack. 

3.  The  stress  at  discontinuity  is  the  long-term  strength. 

4.  In  a controlled  tensile  strain  test,  concrete  is  able  to 
attain  very  large  average  strains. 

5.  The  creep  behavior  of  concrete  in  tension  has  the  same  general 
characteristics  as  in  compression. 

The  concrete  model  used  by  Carino  in  his  experimental  work  was  similar 
to  that  of  Buyukozturk  (ref.  19).  The  method  of  loading  the  specimens 
in  tension  was  similar  to  that  used  by  Kupfer  (ref.  17).  Figure  6 presents 
a plot  of  the  discontinuity  and  failure  envelopes  obtained  from  the  experi- 
ments. Figure  7 shows  the  failure  curves  defined  by  Kupfer  (ref.  17)  and 
Nelissen  (ref.  22).  It  is  interesting  to  note  several  Items  concerning 
Carino's  pair  of  curves.  The  first  is  that  the  discontinuity  curve  is  of 
the  same  shape  as  Kupfer's  ultimate  envelope.  The  second  point  Is  that 
Carino  states  his  ultimate  curve  may  be  a little  Inaccurate  due  to  testing 
apparatus  limitations.  Hence,  he  gives  credence  to  Kupfer's  envelope. 
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The  last  point  to  mention  is  that  under  biaxial  tension  the  discontinuity 
stress  point  and  ultimate  stress  point  essentially  coincide,  whereas 
Kupfer's  criterion  follows  more  along  the  lines  of  brittle  fracture. 
Possibly  the  difference  can  be  accounted  for  by  Kupfer's  using  real  con- 
crete and  Carino  a concrete  model. 


2.  STATICALLY  LOADED  TRIAXIALLY  STRESSED  PLAIN  CONCRETE 


Until  the  last  10  years,  very  little  testing  had  been  done  concerning 
the  triaxial  strengths  of  concrete.  Most  likely  this  was  due  to  the  severe 
requirements  placed  upon  testing  equipment  and  procedures. 

Hannant  (ref.  23)  surveyed  triaxial  experimental  testing  quite  well 
in  his  1968  paper.  From  a plot  of  experimental  data  collected  from  several 
researchers,  he  was  able  to  qualitatively  draw  a failure  envelope  In  tri- 
axial compression.  Figure  8 presents  this  envelope.  The  cross-hatched 
portion  of  this  plot  shows  the  well  defined  biaxial  stress  state. 

Endebrock  and  Traina  of  New  Mexico  State  University  (ref.  24)  concluded 
a biaxial  and  triaxial  testing  program  for  the  Air  Force  in  1972.  Most  of 
the  effort  was  for  triaxial  compression,  with  only  a few  tests  involving 
tension  stresses.  A good  condensation  of  previous  triaxial  tests  up  through 
1969  is  also  presented  in  this  paper.  A general  conclusion  of  this  research 
is  that  the  biaxial  strength  curve  is  quite  well-defined  for  all  principal 
stress  ratios.  Comparisons  of  triaxial  strengths  of  concrete  were  made 
in  pure  compression  only  as  no  previous  tension-compression  results  were 
published.  Figure  9 presents  a typical  plot  of  triaxial  strength  for 
varying  principal  stress  ratios.  Qualitatively  the  plot  confirms  the 
general  shape  of  the  strength  envelope  proposed  by  Hannant  in  Figure  8. 
Further  test  results  by  Terra-Tek  using  cube  specimens  were  incorporated 
into  this  report.  The  results  of  these  tests  provided  strength  curves  of 
similar  shapes;  however,  the  actual  strength  values  were  found  to  be  somewhat 
lower  than  those  of  the  New  Mexico  State  tests.  Figure  10  presents  the 
comparison  of  the  two  tests  for  the  triaxial  compression  tests. 
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EFFECT  OF  DYNAMIC  LOADS  ON  PLAIN  CONCRETE 


The  correct  analysis  for  the  response  of  a complete  structural  sys- 
tem requires  a full  understanding  of  the  material  of  which  the  structure 
is  composed.  This  full  understanding  of  the  material  implies  knowledge 
of  the  material's  simple  behavior  when  experiencing  the  same  type  of  loads 
as  the  structure  of  interest;  i.e.,  rate  and  magnitude  of  loading.  For 
materials  such  as  steel  this  type  of  knowledge  is  well -documented.  Plain 
concrete,  however,  due  to  its  variable  composition,  is  not  well-defined. 
Its  behavior  under  triaxial  quasi-static  stresses  is  just  beginning  to 
become  well -documented.  Also  explanations  or  quantitative  theories  for 
how  concrete  behaves  under  triaxial  quasi-static  stress  vary  from  report 
to  report.  Hence,  until  researchers  begin  to  agree  upon  behavior  patterns 
of  plain  concrete  under  static  stresses,  any  quantitative  or  even  qualita- 
tive explanations  for  the  behavior  of  plain  concrete  under  dynamic  loads 
should  be  scrutlzlned  carefully. 

This  section  presents  a summary  of  work  done  for  determining  behavior 
of  plain  concrete  under  dynamic  loads.  Brink,  et  al  (ref.  25)  conducted 
uniaxial  tests  on  cylindrical  specimens.  The  experimental  study  was  to 
statically  preload  the  plain  concrete  specimens  to  various  stress  levels, 
apply  a small  half-sine  wave  as  a dynamic  pulse,  and  observe  the  behavior. 
The  motivation  for  the  research  was  to  determine  the  extent  that  the 
dynamic  behavior  of  concrete  under  dynamic  loads  would  be  dependent  upon 
the  magnitude  of  preloading.  Three  static  preload  stress  conditions  were 
tested.  The  first  stress,  2000  psi,  represents  the  initial  loading  of 
concrete  where  the  local  stress-strain  are  proportional  to  each  other. 
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The  second  stress,  4000  psi,  corresponds  to  the  nonlinear  portion  of  the 
stress-strain  curve;  and  the  last  preload  stress,  6000  psi,  represented  a 
loading  near  ultimate.  The  lengths  of  the  dynamic  impulses  varied  from 
25  ms  to  2650  ms.  The  results  of  tests  at  the  lowest  preload  stress, 

2000  psi,  showed  that  application  of  dynamic  load  does  not  have  much  effect 
on  the  response  of  concrete.  However,  at  the  higher  preloads,  the  devia- 
tion from  the  normal  stress-strain  curve  was  more  pronounced.  Also,  the 
shorter  the  dynamic  pulse,  the  greater  the  deviation.  A typical  plot  of 
the  results  of  one  of  the  tests  is  given  in  Figure  11.  Figure  12  presents 
a plot  describing  the  increase  in  elastic  modulus  for  each  preload  condition 
as  a function  of  dynamic  pulse  duration.  In  general,  the  study  shows  that 
the  further  the  concrete  is  along  the  stress-strain  curve  the  more  dramatically 
the  stiffness  increases  due  to  dynamic  loads.  Also,  it  was  observed  that 
the  higher  the  loading  rate,  the  more  the  concrete  stiffness  increases. 

Sparks  (ref.  26)  conducted  a study  into  the  effects  of  the  rate  of 
application  of  steadily  increasing  loads  upon  the  static  strength  of  plain 
concrete  in  compression.  Concretes  made  with  gravel,  limestone,  or  lytag 
aggregate  were  tested.  The  results  showed  that  the  static  strength  and 
stiffness  of  concrete  is  increased  at  higher  rates  of  application  of  load. 

It  was  found  that  limestone,  the  stiffest  of  the  aggregates,  produced 
concrete  with  an  improvement  of  only  4%  in  static  strength  with  a hundred- 
fold increase  in  rate  of  loading.  Lytag,  the  least  stiff  of  the  aggregates, 
produced  concrete  with  a strength  improvement  of  16%  for  a similar  increase 
in  rate  of  loading. 


m 


4.  REINFORCEMENT  BOND  SLIP 


Bond  deterioration  between  steel  and  concrete  can  significantly  influence 
the  behavior  of  reinforced  concrete  structures,  especially  under  load  reversals 
(ref.  27). 

Moritz  and  Kaku  (ref.  28)  conducted  an  experimental  program  to  evaluate 
local  bond  slip  deterioration  under  load  reversals.  Figure  13  shows  some 
typical  results  of  their  studies.  They  proposed  a bond  slip  law  but  its 
complexity  makes  it  computationally  undesirable. 
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Figure  7.  Comparison  of  Experimental  Failure  Envelopes  for  Compression- 
Tension  and  Tension-Tension  Loadings 
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Figure  8.  Proposed  Failure  Envelope  for  Concrete 

in  the  Compression  Quadrant  of  Stress  Space 
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Figure  9.  Triaxial  Compression  Strength  Comparisons 


Figure  10.  Normalized  Triaxlal  Strength  Comparison 
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Response  of  Concrete  to  Dynamic  Pulse  with  Preload  Stress 


Static  Prestress 


Rate  of  Loading  Effects  on  the  Increase  In  Stiffness 
and  Strength  of  Concrete 
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SECTION  III 


ANALYTICAL  STUDIES 

In  the  past  decades  there  have  been  many  attempts  to  analyze  and 
comprehend  the  behavior  of  reinforced  concrete  members  and  structural  systems. 
For  various  reasons  a complete  understanding  of  the  subject  has  not  been 
reached.  One  of  the  major  reasons  is  that  the  material  properties  of 
reinforced  concrete  have  not  been  fully  understood.  There  are  insufficient 
data  available  concerning  the  multi  axial  behavior  of  concrete.  Information 
on  the  multiaxial  stress-strain  relationships  is  inadequate.  It  is  also 
not  firmly  established  in  which  condition  the  failure  cracks  initiate  and 
propagate.  Neither  is  the  interaction  of  reinforcing  steel  and  concrete  well 
understood.  Although  there  is  still  a lot  of  uncertainty  regarding  the 
material  property  of  the  reinforced  concrete,  considerable  progress  has  been 
made  recently  in  developing  analytical  models  for  the  reinforced  concrete 
structural  systems,  based  on  the  finite  element  method. 

The  finite  element  method  has  been  proven  as  one  of  the  most  powerful 
and  versatile  tools  in  the  field  of  structural  mechanics.  It  has  been 
described  extensively  in  the  literature.  A comprehensive  discussion  of  the 
theory  and  application  of  the  method  can  be  found  in  any  one  of  the  well-known 
textbooks  (ref.  29).  Only  the  application  of  the  finite  element  method  to 
reinforced  concrete  structures  is  discussed. 

1.  ANALYTICAL  MODELS 

There  are  two  basically  different  approaches  that  are  used  for  the 
application  of  the  finite  element  method  to  reinforced  concrete  structures. 
The  one,  the  discrete  element  approach,  is  to  model  the  concrete  and  the 
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reinforcing  steel  as  discrete  elements.  The  other,  the  composite  element 
approach,  is  to  construct  a composite  modulus  from  individual  properties 
of  concrete  and  reinforcing  steel. 

In  the  discrete  element  approach,  the  reinforced  steel  elements  are 
connected  to  the  concrete  elements  at  the  nodal  points.  Composite  action 
is  enforced  through  compatibility  at  the  nodes.  A bond  linkage  element 
may  be  used  to  simulate  bond  slip  behavior. 

Ngo  and  Scordelis  (ref.  30)  presented  the  first  finite  element  model 
for  the  analysis  of  a reinforced  concrete  beam  with  the  discrete  element 
approach.  Simple  beams  were  represented  by  two-dimensional  triangular 
finite  elements. 


Ngo,  Scordelis  and  Franklin  (ref.  31)  extended  the  work  by  using  refined 
linkage  elements  to  include  the  influence  of  stirrups,  dowel  shear,  aggre- 
gate Interlock,  etc.  Nilson  (ref.  32)  also  refined  the  work  of  Ngo  and 
Scordelis  by  introducing  nonlinear  material  properties  and  nonlinear  bond- 
slip  relationships  into  the  model  to  analyze  beam-type  structures.  Taylor, 
etc.  (ref.  33,  34)  adopted  Kupfer,  Hilsdorf  and  Rusch  (ref.  35)  experimen- 
tal results  incorporated  into  the  discrete  element  approach  to  predict  the 
behavior  of  axisynmetrlc-type  reinforced  and  prestressed  concrete  structures. 
Salam  & Mohraz  (ref.  36)  applied  the  discrete  element  approach  to  study 
the  nonlinear  behavior  of  planar  reinforced  concrete  structures  with  octa- 
hedral shear  stress  failure  criteria  to  check  crack  or  yield.  Since  1973, 
Argyris  et  al  (refs.  37,  38,  39)  studied  nonlinear  behavior  and  ultimate 
load  for  prestressed  and  reinforced  concrete  reactor  structures  by  the  dis- 
crete element  approach.  Litton  (ref.  27)  developed  a two-dimensional  finite 
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element  for  predicting  brittle  material  response  under  cyclic  loading, 
which  included  crack  formation,  closure,  and  reopening  effects.  In  addi- 
tion, a bond  link  element  models  concrete- reinforcement  bond  and  slip. 

The  results  compared  favorably  with  experimental  data. 

The  more  recent  work  by  Same  (ref.  40)  uses  three-dimensional  finite 
element  analysis  for  reinforced  and  prestressed  concrete  structures  by 
using  the  discrete  element  approach.  For  the  discrete  element  approach 
applied  to  reinforced  concrete  structures,  concrete  is  idealized  as  2-D 
or  3-D  elements,  and  the  reinforcing  steel  is  modeled  by  2-D  elements  or 
truss-bar  elements,  as  shown  in  Figure  14. 


The  composite  element  approach  combines  the  variation  in  steel  and  concrete 
properties  and  constructs  a composite  modulus  or  stiffness  for  a reinforced 
concrete  element  to  account  for  the  composite  behavior  of  the  element. 

Usually  the  bond  between  reinforcing  steel  and  concrete  is  assumed  to  be 
fully  effective.  But  some  Investigators  (ref.  41)  still  included  the 
influence  of  the  bond  slip  in  the  construction  of  the  composite  modulus. 


Mikkola  and  Schnobrich  (ref.  42)  applied  the  finite  element  method  to 
study  the  response  of  reinforced  concrete  shells  by  constructing  a composite 
modulus  as  follows: 
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ac  • -£c  + as 


(4) 


in  which  C 


material  property  matrix  for  a composite  element 
material  property  matrix  for  reinforced  steel 
material  property  matrix  for  plain  concrete 
relative  area  of  the  reinforcing  steel 
A$/A 
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A$  = steel  area 

A = total  area 

ac  = relative  area  of  the  plain  concrete 

= 1 - a$ 

Scanlon  (ref.  43)  analyzed  slabs  with  the  composite  element  approach. 
Time  dependent  effects  of  shrinkage  and  creep  were  also  included. 

Yuzugullu  and  Schnobrich  (refs.  44,  45)  obtained  the  composite  material 
property  matrix  for  the  reinforced  concrete  model  as  the  sum  of  the  material 
property  matrices  of  concrete  and  reinforcement: 

C - Cc  . l «,  cj  <W 

in  which  o^.  = the  percentage  of  reinforcement  in  i -direction 

_ area  of  reinforcement  in  i -direction 
concrete  area 

Any  number  of  reinforcement  directions  can  be  accommodated.  Due  to 
the  presence  of  reinforcement,  the  reinforced  concrete  model  is  no  longer 
isotropic.  The  directions  of  the  principal  stresses  and  the  principal 
strains  would  be  different.  Three  types  of  structures,  a deep  beam,  a 
shear  panel,  and  a shear  wall-frame,  were  analyzed.  Aggregate  interlock, 
or  the  ability  of  a cracked  element  to  transfer  shear,  was  included  In  the 
analysis.  A constant  value  of  shear  stiffness  was  assumed  for  all  the 
cracked  elements. 
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Cervenka  (ref.  46),  Cervenka  and  Gerstle  (refs.  47,  43)  studied  the 
in-plane  response  of  a wall  panel  with  the  composite  material  property 
matrix  calculated  in  the  same  way.  The  shear  resistance  of  the  cracked 
concrete  owing  to  the  aggregate  interlocking  and  reinforcement  dowel 
effect  was  neglected.  This  may  contribute  to  the  underestimation  of 
the  stiffness  of  the  theoretical  model  in  the  cracked  stage. 

Hand,  Pecknold  and  Schnobrich  (refs.  49,  50)  studied  the  load-deflection 
history  up  to  failure  of  reinforced  concrete  slabs  and  shells  by  a layered 
system  (Figure  15).  The  element  is  divided  into  a number  of  layers  over 
the  depth.  Each  layer  is  assumed  to  be  in  a state  of  plane  stress  and 
may  have  different  material  properties  corresponding  to  its  material  state. 

Lin  (ref.  51),  Lin  and  Scordelis  (refs.  52,  53)  used  the  same  approach  to 
do  nonlinear  analysis  of  reinforced  concrete  shells  of  general  form.  The 
differences  between  the  two  works  were 

1.  A flat  triangular  element  adaptable  to  free-form  shell  surfaces  and 
arbitrary  boundary  geometries  was  used  in  Lin's  work,  while  Hand's 
work  used  a shallow  shell  element  to  analyze  plates  and  shallow 
shells. 

2.  The  tension  stiffening  effect  of  the  concrete  between  cracks, 
neglected  In  Hand's  work,  was  Included  in  Lin's  work.  It  was 
found  in  Lin's  work  that  the  effect  had  significant  influence  on 

§ 

the  post-cracking  load-deflection  response  of  underreinforced 
concrete  structures. 


Wanchoo  and  May  (ref.  54)  also  adopted  a similar  approach  to  performing 
cracking  analysis  of  reinforced  concrete  plates.  It  may  be  noted  that 
the  classical  Kirchhoff  assumption  was  incorporated  into  the  layer  system 


approach.  The  assumption  states  that  a normal  remains  normal  during  deforma- 
tion, so  each  layer  Is  considered  to  be  in  a state  of  plane  stress  with  prop- 
erties defined  through  biaxial  stress-strain  relationships.  Because  of 
this  assumption,  shear  modes  of  failure  with  their  associated  diagonal  crack- 
ing cannot  be  solved  with  this  scheme.  Three-dimensional  models  become  a 
necessity  for  the  investigation  of  any  failure  development  where  shear  is  a 
major  consideration. 

Rajagopal  (ref.  55)  adopted  the  same  approach  to  study  the  reinforced 
concrete  beams,  beam  columns  and  slabs.  Adlstedt  (ref.  56)  used  a similar 
approach  to  investigate  the  nonlinear  behavior  of  reinforced  concrete 
frames.  The  effects  of  geometric  nonlinearities  were  considered  in  both 
works. 

Adham,  Braunick  and  Isenberg  (ref.  41)  developed  for  AFWL  an  equivalent 
continuum  approach,  similar  to  a composite  element  approach  applicable  to 
beams,  beam  columns  and  ring  beams  subjected  to  static  and  dynamic  loading. 

A single  stress-strain  relationship  was  derived  directly  for  reinforced 
concrete.  Bond  slip  and  dowel  shear  relations  were  incorporated  into  the 
model . 

Basically  speaking,  the  discrete  element  approach  is  more  straightforward 
in  the  mathematical  procedures  and  simpler  for  computer  programming.  The 
composite  element  approach  is  more  flexible  for  different  arrangements  of 
the  reinforcement.  The  orientation  and  location  of  the  reinforcement  do 
not  influence  the  arrangement  of  element  mesh  for  the  composite  element 
approach  so  much  as  it  does  for  the  discrete  element  approach. 


So  far,  the  finite  element  method  has  provided  the  basis  for  all  the 
analytical  solutions,  and  the  nonlinear  effect  is  based  on  the  basic  non- 
linear stress-strain  law.  An  alternative  approach,  the  modified  El  approach, 
may  replace  the  composite  element  approach  to  study  the  nonlinear  response  of 
reinforced  concrete  structures.  The  approach  is  based  on  an  empirical 
moment-curvature  relationship,  wherein  different  flexural  rigidities  are 
assumed  for  different  material  states.  Jofrlet  and  McNiece  (ref.  57)  used 
a bilinear  moment-curvature  relation  to  study  the  reinforced  concrete  plate 
bending  problem.  Bell  (ref.  58)  adopted  a similar  material  idealization  to 
study  slabs  and  shells.  Cracking  of  flat  slabs  was  studied  with  this  approach 
by  Jofriet  (ref.  59).  The  advantage  in  using  the  modified  El  approach  would 
be  in  the  reduction  of  computer  costs  for  similar  analyses. 

2.  MATHEMATICAL  MATERIAL  MODELS 

No  matter  what  analytical  approach  is  used,  it  requires  a mathematical 
material  model  for  each  individual  component,  concrete  and  reinforcing 
steel  to  simulate  the  material  behavior.  Compared  to  concrete,  the  material 
property  of  steel  is  very  well-defined.  Current  yield  criterion  and  stress- 
strain  relations  for  steel  have  been  demonstrated  to  be  adequate  in  predicting 
steel  behavior  even  in  the  plastic  range.  In  contrast,  the  behavior  of  con- 
crete is  complicated » and  there  still  exists  a lot  of  uncertainty  regarding 
the  material  property  of  concrete.  Any  mathematical  material  model  applied 
in  analytical  approaches  to  simulate  a realistic  behavior  of  concrete  must 
employ  drastic  idealizations.  Many  models  have  been  proposed  to  simulate 
reinforced  concrete  behavior  using  finite  element  methods.  Review  of  these 
models  is  discussed  in  the  following  paragraphs. 
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a.  Linear  Elastic  Model  with  Linkage  Element 

The  first  finite  element  model  for  the  analysis  of  a reinforced 
concrete  beam  was  developed  by  Ngo  and  Scordells  (ref.  30).  In  their 
studies,  simple  beams  were  analyzed  in  which  the  concrete  and  steel 
reinforcement  were  represented  by  plane  stress  triangular  finite 
elements.  Special  linkage  elements  were  developed  to  connect  the 
reinforcing  steel  to  the  concrete  to  represent  the  bond-slip  phenomena 
(Figure  16).  Linear  elastic  analyses  were  performed  on  reinforced 
concrete  beams  with  predicted  crack  patterns  to  determine  principal 
stresses  in  the  concrete,  stresses  in  the  steel  reinforcement,  and  bond 
stresses.  Ngo,  Scordelis,  and  Franklin  (ref.  31)  extended  this  work 
to  study  shear  in  beams  with  diagonal  tension  cracks.  The  model  was 
refined  by  representing  the  steel  reinforcement  and  the  concrete  by 
quadrilateral  elements,  and  the  vertical  stirrups  were  modeled  by  truss 
bar  elements  (Figure  17).  The  influence  of  dowel  shear,  aggregate 
interlock,  and  horizontal  splitting  along  reinforcement  near  the  support 
was  considered. 

However,  a linear  elastic  model  is  not  suitable  for  the  study  of 
the  behavior  of  reinforced  concrete  structures  over  the  entire  range 
of  loading.  Nilson  (ref.  32)  introduced  nonlinear  material  properties 
and  a nonlinear  bond-slip  relationship  into  the  model  and  used  an 
incremental  loading  technique  to  account  for  these  nonlinearities. 

In  all  these  studies,  cracking  of  concrete  was  Introduced  by  physically 
splitting  a node  and  allowing  a separation  to  develop  along  interelement 
boundaries.  Cracks  were  therefore  restricted  to  propagate  only  along 
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Interelement  boundaries.  In  order  to  be  effective  and  have  the  element 
boundaries  in  the  proper  location  for  the  crack,  previous  knowledge  of 
the  location  of  the  crack  was  required,  or  undue  constraint  was  placed  on 
the  direction  of  the  development  of  the  crack  path  (Fiqure  18).  Furthermore, 
the  disruption  of  structural  continuity  due  to  cracks  results  In  a 
continual  change  in  the  topology  of  the  structure.  This  presents  a number 
of  analytical  difficulties  and  some  computational  inefficiency  and 
restricts  the  model  application  to  static  load  cases  only. 

Recently,  Ngo  (ref.  60)  developed  a network- topological  approach  to 
simulate  crack  growth  by  automatically  generating  crack  line  In  the  finite 
element  model.  However,  this  approach  is  still  in  the  infant  stage,  and 
research  is  needed  to  make  this  model  more  workable. 

b.  Orthotropic  Material  Model 

Rashid  (ref.  61)  performed  a nonlinear  analysis  of  axisymmetric 
prestressed  concrete  nuclear  reactor  pressure  vessels  with  an  ortho- 
tropic material  model.  Instead  of  redefining  the  new  topology  of  the 
structure  (i.e.,  new  finite  element  mesh  and  connectivities),  the 
cracking  of  concrete  was  taken  Into  account  by  modifying  the  material 
properties  and  redistributing  the  unbalanced  (residual)  stresses.  The 
mechanics  of  accounting  for  tension  cracking  Involved  computation  of 
principal  stresses  or  strains  which  were  then  compared  with  a cracking 
criterion.  If  the  criterion  was  met  or  exceeded,  a crack  was  assumed 
to  have  developed,  stiffness/modulus  normal  to  the  crack  was  reduced 
to  zero,  and  the  stress  field  associated  with  the  cracking  was  redis- 
tributed to  the  rest  of  the  structure.  This  model  allows  free  development 
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of  cracking  without  predefining  the  cracking  path,  and  permits  the 
use  of  the  same  structural  topology  throughout  the  solution,  rather 
than  the  necessity  for  splitting  the  nodes  after  cracking  and 
establishing  a new  topology.  However,  none  of  the  nonlinear  effects 
of  concrete  in  the  compression  region  were  considered  realistically. 

c.  Hypo-Elastic  Model  (Nonlinear  Elastic  Model) 

Liu,  Nilson  and  Slate  (ref.  62)  developed  an  equivalent  constitu- 
tive relationship  for  concrete  based  on  an  orthotropic  material  assump- 
tion using  some  experimental  data  of  biaxial  compression  tests.  The 
equivalent  stress-strain  relationship  is  given  as 
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stress  and  strain  respectively 
Poisson's  ratio 

ratio  of  the  principal  stress  In  the  orthogonal 
direction  to  the  principal  stress  in  the  direction 
considered. 

Initial  tangent  modulus  in  uniaxial  loading 
current  modulus  at  the  ultimate  load 


For  concrete  in  tension.  It  was  assumed  that  the  behavior  is  linear 
up  to  failure.  For  an  elastic  orthotropic  material  In  biaxial  stress, 
the  constitutive  relationship  may  be  written  as 


39 


r 


* 


“ — 

°1 

Ei*  v* 

el 

< °2 

> = 

V*  2* 

-< 

e2 

> (7) 

T12 

G * 

1 2 

Yl2 

- _ 

L J 

where 


E * 
1 


Ei 

1-^  E2/Ex 


E * = * 

l-vf  Ej/Ej 


G * = 
U1  2 


Ei  E j 


Ei+Ez^VjE- 


in  which 


ai»  a2  » Gl»  £2 

T12»  Y12 


E 

1 » 


E 

2 


principal  stresses  and  strains,  respectively, 
shear  stress  and  strain  in  principal  direction, 
respectively. 

uniaxial  tangent  modulus  in  principal  direction  1 
and  2,  respectively. 

uniaxial  Poisson's  ratio  in  the  principal 
direction  1. 


Kupfer  and  Gerstle  (ref.  63)  developed  another  hypo-elastic  model 
based  on  an  isotropic  formulation  using  variable  bulk  moduli  and  shear 
moduli.  Another  hypo-elastic  model  for  triaxial  states  of  loading  was 
derived  by  Coon  and  Evans  (ref.  64).  Incremental  constitutive  laws 
and  their  associated  failure  criteria  were  presented  for  plain  concrete 
based  on  an  assumption  of  isotropic  material.  Argyris,  et  al. 

(ref.  37)  discuss  some  of  the  problems  Involved  in  using  the  hypo- 
elastic  model.  With  this  model  the  modulus  E becomes  anisotropic 
in  the  nonlinear  range,  even  though  the  behavior  Is  isotropic  in  the 
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Initial  linear  range.  This  Implies  different  behavior  in  the  principal 
stress  directions.  It  also  suggests  that  the  principal  directions  of 
stress  and  strain  are  different.  Introducing  coupling  between  normal 
stresses  and  shear  strains.  It  requires  an  Increase  in  the  number 
of  material  moduli  to  define  this  anisotropic  behavior.  In  addition, 
failure  condition  must  be  supplied  to  match  the  realistic  deformation 
behavior  of  concrete  in  the  ultimate  response  regime. 

d.  Combined  Behavior  Model 

In  the  last  several  years  research  has  been  directed  primarily 
toward  development  of  more  realistic  analytical  models  for  reinforced 
concrete  and/or  plain  concrete.  Most  of  them  may  be  included  in  this 
group. 

A typical  stress-strain  relation  for  plain  concrete  under  uniaxial 
condition  is  shown  in  Figure  19.  The  linear  elastic  region  is  followed 
by  a hardening  range  in  which  Internal  microcracking  progressively 
weakens  the  material  property  up  to  the  ultimate  stress  point.  Upon 
continued  loading  the  curve  deteriorates,  indicating  softening  in  the 
structure  up  to  the  rupture  point  at  which  sudden  collapse  is  observed. 
The  behavior  in  tension  is  linear  elastic  followed  by  a sudden  crack. 

Cervenka  (refs.  46,  47,  48)  studied  inelastic  behavior  of  reinforced 
concrete  panels  by  assuming  reinforcing  steel  to  be  elastic  perfectly 
plastic  in  both  compression  and  tension  (Figure  20).  The  concrete  was 
assumed  to  be  elastic  perfectly  plastic  in  compression  and  elastic 
brittle  in  tension  (Figure  21).  The  von  Mises  yield  criterion  for 
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compression  plasticity  and  maximum  normal  stress  theory  for  tension 
cracking  was  adopted  to  approximate  the  failure  surface  of  the  con- 
crete in  biaxial  stress  state  (Figure  22).  Yuzugullu  and  Schnobrich 
(refs.  44,  45)  applied  the  same  material  model  for  the  prediction  of 
inelastic  behavior  of  shear  wall -frame  systems. 

Hand,  Pecknold  and  Schnobrich  (refs.  49,  50)  derived  a layered 
system  to  analyze  reinforced  concrete  plates  and  shallow  shells.  The 
concrete  was  considered  as  bilinear  elastic  or  perfectly  plastic  material 
in  compression  and  elastic  brittle  in  tension  (Figure  23).  A yield 
criterion  in  biaxial  stress  state  proposed  by  Kupfer,  Hllsdorf  and 
Rusch  (ref.  35)  was  used.  The  yield  surface  was  modified  to  serve  as 
the  transition  criterion  to  determine  the  boundaries  of  the  different 
material  property  zones.  This  was  accomplished  by  simply  scaling  down 
the  yield  surface  of  Kupfer,  Hilsdorf  and  Rusch.  A crush  surface, 
analogous  to  the  yield  surface  but  In  terms  of  strain,  was  proposed 
to  limit  the  deformation  behavior  of  the  compression  concrete  (Figure  24). 

Lin  and  Scordelis  (refs.  51,  52,  53)  proceeded  along  a similar 
approach  to  perform  nonlinear  anlaysls  of  reinforced  concrete  shells  of 
general  form.  However,  the  material  model  and  failure  conditions  were 
slightly  different  from  Hand's  work.  The  concrete  was  assumed  to  be 
elastic  perfectly  plastic  in  compression  and  elastic  brittle  in  tension 
(Figure  25).  Because  of  the  tension  stiffening  effect,  the  concrete 
was  considered  to  release  its  stress  gradually  after  it  cracked.  There- 
fore, there  existed  a cracked  unloaded  portion  in  the  tension  region  of 
the  stress-strain  curve  for  concrete  (Figure  26).  The  von  Mises  yield 
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criterion  was  used  to  approximate  the  failure  surface  in  biaxial  com- 
pression (Figure  27).  The  associated  flow  rule  was  assumed  to  govern 
the  post-yielding  stress-strain  relations  for  concrete.  A crush  sur- 
face, analogous  to  the  yield  surface  but  in  terms  of  strains,  was 
postulated  to  define  the  complete  crush  for  the  yielded  concrete 
(Figure  28).  The  maximum  stress  failure  theory  was  adopted  for  biaxial 
tensile  stress  state.  For  the  tension-compression  stress  state,  a straight 
line  was  used  to  approximate  the  true  failure  surface  (Figure  27).  The 
reinforcing  steel  in  both  Lin’s  and  Hand's  works  was  considered  as  an 
elastic  perfectly  plastic  material.  In  addition,  the  reinforcing  bars 
were  assumed  to  carry  only  axial  stress. 

Gormack  (ref.  65)  adopted  a material  model  similar  to  the  one  proposed 
by  Lin  and  Scordelis  to  analyze  shell  walls  and  two-dimensional  structures. 
The  only  difference  was  that  the  stress-strain  curve  for  concrete  in 
compression  proposed  by  Liu,  Nilson  and  Slate  (ref.  62)  was  used  instead 
of  the  assumption  of  elastic  perfectly  plastic  behavior. 

Salem  and  Mohraz  (ref.  36)  analyzed  planar  reinforced  concrete 
structures  with  finite  elements.  The  stress  strain  curve  for  concrete 
with  an  artificial  unloading  portion  in  tension  and  elastic  perfectly 
plastic  behavior  in  compression  was  assumed  (Figure  29).  The  octahedral 
shear  stress  failure  criterion  was  considered  to  check  yielding  or 
cracking  (Figure  30).  Crushing  of  concrete  took  place  If  the  equivalent 
uniaxial  plastic  strain  exceeded  a prescribed  value.  The  reinforcement 
was  treated  as  an  elastic  perfectly  plastic  material  with  von  Mlses 
yielding  criterion. 


Darwin  and  Pecknold  (ref.  66)  Introduced  an  analytical  model  for  cyclic 


biaxial  loading  for  reinforced  concrete.  The  constitutive  model  for 
plain  concrete  was  developed  to  match  experimental  data  on  monotonic 
biaxial  and  cyclic  uniaxial  loading,  and  extended  to  cover  cyclic 
biaxial  loading  and  multidirectional  cracking.  The  analytical  biaxial 
strength  envelope  for  the  proposed  model  is  snown  in  Figure  31.  The 
proposed  model  under  cyclic  load  is  given  in  Figure  32.  The  stress- 
strain  equation  suggested  by  Saenz  (ref.  67,  68)  was  used  for  the 
ascending  portion  of  the  stress-strain  curve  for  the  proposed  model 
(Figure  33).  An  envelope  curve,  a set  of  coirmon  points,  and  a set  of 
turning  points  were  Introduced  to  define  unloading  and  reloading  behaviors 
of  concrete  under  cyclic  loading.  These  introductions  made  the  proposed 
model  more  complicated  in  analytical  procedure. 

Same  (ref.  40)  extended  the  analytical  method  to  perform  three- 
dimensional  nonlinear  analysis  of  reinforced  and  prestressed  concrete 
structures.  The  uniaxial  stress-strain  relation  for  concrete  was  based 
on  the  one  proposed  by  Saenz  (Figure  34).  Concrete  was  checked  for 
tensile  failure  using  the  stress  in  the  principal  direction.  The  maxi- 
mum stress  il are  theory  was  adopted  for  triaxlal  tensile  stress  state. 
The  tensil.  strength  was  assumed  to  decrease  linearly  with  increasing 
compressive  stress  in  other  principal  directions.  The  tensile  failure 
surface  is  shown  in  Figure  35.  Compressive  failure  of  concrete  was  also 
checked  using  principal  stresses.  The  ultimate  strength  of  concrete 
was  assumed  to  be  affected  by  compressive  stresses  In  the  other  principal 
directions.  The  failure  surface  proposed  by  Khan  and  Sangy  (ref.  69) 
was  used  because  of  better  compatibility  with  biaxial  failure  criterion. 


The  surface  was  a cone  center  along  the  average  axis  of  all  three 
principal  compressive  stress  directions  (Figure  36).  No  flow  rate 
was  used  for  the  plasticity  of  concrete  in  compression. 
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The  modulus  of  elasticity  of  concrete  was  determined  based  on  the 
I 

uniaxial  stress-strain  relationship.  Each  principal  stress  was  taken 
independently  to  determine  the  reduction  in  the  modulus  in  that  direc- 
tion due  to  the  effect  of  the  stresses  in  the  other  principal  directions. 

Unloading  of  concrete  in  compression  was  treated  differently  than 
that  in  tension.  The  initial  modulus  was  used  for  unloading  from  a state 
of  compression.  During  reloading,  the  initial  modulus  was  used  until 
the  state  when  the  compressive  stress  in  the  concrete  reached  its 
previous  value.  In  tension  while  unloading,  if  the  total  strain  across 
the  crack  was  still  positive,  the  crack  was  considered  to  be  open  and 
no  rigidity  was  considered  for  that  direction.  If  the  strain  now  was 
negative,  the  crack  was  closed  and  the  rigidity  of  concrete  in  compression 
was  used,  but  this  principal  direction  was  fixed.  For  reloading  when 
tension  strain  was  developed,  the  same  crack  would  open  again. 

The  reinforcing  steel  was  modeled  as  an  elasto-plastic  material 
having  the  yield  and  ultimate  properties  controlled  by  von  Mises 
criterion.  The  recent  work  done  by  Rajagopal  (ref.  55)  may  be  included 
in  this  combined  behavior  model. 

Uncracked  concrete  element  with  this  type  of  modeled  behavior  was 
considered  as  an  Isotropic  material.  After  cracking,  the  element  was 
considered  as  an  anisotropic  material  with  nonlinear  response  properties. 
Stresses  and/or  strains  in  the  principal  directions  were  used  to 
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check  the  nonlinear  behavior  of  concrete  In  compression.  The  consti- 
tutive matrix  was  expressed  In  the  local  coordinate  system  of  principal 
stresses  and  strains.  Once  the  material  was  considered  anisotropic, 
the  directions  of  the  principal  stresses  and  strains  were  no  longer 
coincident  with  the  coordinate  axes.  This  condition  resulted  in  the 
coupling  effect  between  normal  stresses  and  shea'-*  strains.  This 
coupling  effect  was  always  neglected  by  most  investigators  using  the 
combined  behavior  model.  Since  the  principal  stresses  were  used  to 
check  the  formation  of  cracks,  the  directions  of  cracks  were  orthogonal 
to  each  other.  This  assumption  is  valid  if  cracks  occur  at  the  same 
time.  Cracks  developing  at  different  loading  stages  are  not  necessarily 
orthogonal  because  of  the  change  of  state  of  stress  in  the  section. 

This  type  of  behavior  was  also  not  considered  by  most  researchers 
using  the  combined  behavior  model.  Few  have  made  some  modifications  by 
assuming  that  a second  crack  develops  in  the  direction  with  an  oblique 
angle  to  the  first  crack,  but  it  doesn't  necessarily  match  the  actual 
situation.  In  addition,  for  cyclic  loading  cases,  it  becomes  necessary 
to  keep  track  of  the  direction  and  width  of  each  crack  in  the  structure, 
which  is  computationally  complex  and  expensive. 

Elastic  Brittle  Model 

Litton  (ref.  27)  adopted  an  elastic  brittle  material  model  for  plain 
concrete  to  study  a reinforced  concrete  beam  and  a wall  panel  under 
cyclic  loading.  The  emphasis  was  placed  on  tensile  cracking,  with  no 
attention  paid  to  the  nonlinear  compression  response.  Hence  a simpli- 
fied fracture  surface  with  a constant  fracture  stress  was  assumed.  A 


concise  and  general  constitutive  theory  was  presented  which  permits 
continual  monitoring  of  crack  behavior,  including  crack  closure  and 
reopening  under  cyclic  load.  A bond  link  element  was  also  developed. 

f.  Linear  Elastic  Plastic  Strain  Hardening  and  Fracture  Model 

Chen  and  Chen  (refs.  70  and  71)  studied  some  punch  indentation  problems 
of  concrete  by  proposing  constitutive  relations  for  behavior  of  concrete. 
The  concrete  was  considered  to  be  a linear  elastic  plastic  strain- 
hardening and  fracture  material.  The  nonlinear  effects  caused  by  the 
plasticity  of  concrete  in  compression  and  cracking  of  the  concrete 
in  tension  were  considered  in  the  formulation  of  these  equations.  The 
constitutive  model  was  based  on  three  fundamental  assumptions;  the  shape 
of  an  initial  discontinuous  surface,  the  evolution  of  subsequent  loading 
surfaces,  and  the  formulation  of  an  appropriate  hardening  rule.  In  other 
words,  initial  discontinuous  surface,  failure  surface,  and  stress-strain 
curve  (Figure  37)  shall  be  defined  before  the  constitutive  model  is 
constructed.  The  failure  surface  was  defined  in  the  stress  space  such 
that  once  the  stress  state  reached  this  surface,  the  material  would 
completely  rupture  and  could  not  resist  any  further  loading.  The  initial 
discontinuous  surface  was  a stress  surface  that  could  be  reached  only 
by  elastic  action.  If  concrete  was  stressed  beyond  its  initial  discon- 
tinuous surface,  microcracks  would  propagate  in  the  concrete  resulting 
in  irrecoverable  deformations  upon  unloading.  If  concrete  was  reloaded 
after  such  an  unloading  had  occurred,  additional  Irrecoverable  deforma- 
tions would  not  occur  until  this  new  subsequent  discontinuous  surface  was 
reached. 
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Concrete  failure  criterion  depends  on  deviatoric  stresses  and  hydro- 
static pressure.  Two  different  but  similar  functions  were  used  for  the 
failure  surface  In  the  compression  region  and  in  the  tension-compression 
region;  i.e.. 
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where  the  positive  sign  in  the  third  term  defines  the  compression  region, 
and  negative  sign  defines  the  tension-compression  region.  The  failure 
criterion  is  a function  of  Ix , the  first  invariant  of  the  stress  state 
tensor.  The  dependence  of  failure  criterion  on  J3,  the  third  invariant 
of  the  deviatoric  stress  tensor  , was  not  considered  to  ensure  the 
simplicity  of  the  formulation  of  the  constitutive  relations.  The 
inclusion  of  J3  will  definitely  improve  the  correlations  between  the 
model  and  experimental  results. 

The  initial  discontinuous  surface  Is  assumed  to  take  the  same  form 
as  that  of  the  failure  surface. 


fo^ai j ^ = 3 J2  ' 36  l'2  f T2  V + 3 Vi  = Tc 


(9) 


in  which  Aq,  tq,  A^  and  x^  are  material  constants,  which  can  be  determined 
from  simple  tests. 

The  subsequent  loading  surfaces  are  assumed  to  be  similar  in  form  to 
those  of  the  initial  discontinuous  and  failure  surfaces. 
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If  plotted  in  principal  stress  space,  the  Initial  discontinuous 
and  failure  surfaces  are  given  (Figure  38).  Good  agreement  with  experi- 
mental results  in  the  biaxial  cases  [i.e.,  Kupfer,  etc.  (ref.  35)  and 
Vllle  (ref.  72)]  are  also  shown  in  Figure  39. 

The  elastic-plastic  stress-strain  incremental  relationships  of 
concrete  are  derived  using  the  classical  theory  of  plasticity.  The 
general  form  of  the  incremental  elastic-plastic  strain-stress  relation- 
ship may  be  obtained  in  the  form 


C1jkl  dekl 


(10) 


The  stress-strain  relations  expressed  explicitly  in  matrix  form  may 
be  used  in  finite  element  analysis. 

The  criterion  for  unloading  from  plastic  state  is  given  as 


df 


3f 


to 


1j 


dou  < 0 


and 

f(CTij>  = T 

The  criterion  for  loading  is 


OD 


(12) 


df 


3f 

3oij 


(13) 


or  df 


df 

301j 


> 0 


(14) 
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It  may  be  noted  that  no  consideration  would  be  given  by  this  model 


to  determine  the  orientation  and  the  propagation  of  each  individual 
crack.  The  microcracks  were  assumed  to  be  randomly  oriented.  The 
cracking  was  considered  from  a macroscopic  point  of  view  which  could 
be  simulated  by  the  behavior  of  an  equivalent  continuum.  The  change 
of  the  material  properties  due  to  cracking  could  be  considered  by 
changing  constitutive  relations  at  integral  points. 

In  addition,  the  Incremental  stress-strain  relations  were  derived 
based  on  the  classical  elastic-plastic  model.  It  has  been  in  this 
field  that  the  finite  element  method  has  been  applied  most  successfully. 

g.  Other  Material  Models 

Wiliam  and  Warnke  (ref.  73)  developed  a three-parameter  failure 
surface  for  concrete  subjected  to  trlaxial  loading  in  the  tensile  region 
and  in  the  region  of  low  compressive  stresses  (Figure  40).  The  failure 
surface  was  convex,  continuous  with  a tension  cut-off  in  the  tension 
regime.  The  model  was  refined  by  adding  two  additional  parameters  for 
describing  curved  meridians,  thus  extending  the  range  of  application  to 
the  high  compression  region  (Figure  41).  The  constitutive  relations 
were  constructed  based  on  an  elastic  plastic  formulation  with  a brittle 
failure  condition  in  the  tensile  regime.  Based  on  the  normality  principal, 
explicit  expressions  were  developed  for  the  inelastic  deformation  rate 
and  the  corresponding  incremental  stress-strain  relation.  No  consideration 
was  given  to  the  formation  of  each  individual  crack.  Cracking  was  con- 
sidered from  a macroscopic  point  of  view  which  could  be  simulated 
by  the  behavior  of  an  equivalent  continuum.  It  Is  a concept  similar  to 
the  one  proposed  by  Chen  and  Chen. 


Green  and  Swanson  (ref.  74)  derived  a three-dimensional  static 
constitutive  relation  for  concrete  based  on  an  Initial  plasticity 
model  used  by  Swanson  for  rock,  with  the  yield  cap  suggested  by  Baron, 
et  al.  An  associated  flow  law  and  a segmented  yield  function  have  been 
used  with  a hardening  law  which  Includes  the  shear  and  dilatation  terms. 
But  it  cannot  handle  unloading,  nor  does  it  fit  in  detail  the  stress 
dependence  of  yield.  It  also  cannot  uniquely  represent  the  maximum 
stress  surface.  Like  other  plasticity  cap  models,  it  appears  too  compli- 
cated to  fit  data  to  it. 

Bazant  (ref.  75)  extended  the  endochronic  theory  for  metals  (ref.  76) 
to  describe  the  inelastic  behavior  and  failure  of  concrete.  The  model 
allows  a far  better  description  of  the  behavior  of  concrete  than  other 
constitutive  models.  It  gives  an  apparently  complete  description  of 
the  nonlinear  behavior  of  concrete,  including  creep.  Although  the 
theory  suitably  fits  finite  element  analysis,  incorporation  would 
require  considerable  modification  to  current  programs  available. 
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Figure  14.  Discrete  Element  Model 


Figure  15.  Composite  Element  Model 
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Figure  19.  Stress-Strain  Relationship. 

Uniaxial  Behavior  of  Plain  Concrete 
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tension 

compression 


unloading 


Figure  20.  Uniaxial  Stress-Strain  Relation  for  Steel  Reinforcement 


compression 


unloading 


Figure  21.  Uniaxial  Stress-Strain  Relation  for  Concrete 


Figure  22.  Biaxial  Strength  of  Concrete 
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Figure  24.  Boundaries  between  the  Different  Material  Behavior  Real 
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Figure  25.  Idealized  Uniax 
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Figure  26.  Gradual  U 
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Figure  29.  Actual  and  Assumed  Uniaxial  Stress-Strain  Curve  for  Concrete 


I 


a 


Where:  f,e  Concrete  stress  and  strain,  respectively 

eQ  Ultimate  strain  = (31.5  - Vf£)*  VT£  * 10“ 3 
E Initial  Young's  modulus 

Eq  Secant  modulus  at  the  ultimate  stress  (f^/e  ) 


Figure  34.  Uniaxial  Stress-Strain  Relationship  for  Concrete 


+a 


where: 


n = the  ratio  of  the  biaxial 
compressive  strength  to  the 
uniaxial  strength 
= 1.25 

11  = ai  + °2  + a3 

12  = f [(a1-a)2+(a2-5)2+(a3-a)2] 


and  a = 


Oj+  o2+  a3 


'in  which  aj,a2,a3  are  principal 
stresses. 


Figure  36.  Triaxial  Compressive  Stresses  Failure  Surface 
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Figure  37.  Typical  Stress-Strain  Curve  for  Concrete 
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Figure  39. 


Failure  and  Initial  Discontinuous  Curves  In  Biaxial 
Principal  Stress  Space:  (a)  Experimental  Data  from 

’>);  (b) 


Kupfer  et  al . (ref.  35) : 
Ville  (ref. 


Experimental  Data  from 


Strength  Ratios  « = 1.3,  a = 0.10 


SECTION  IV 


CONSTITUTIVE  MODELS 

Based  on  the  study  presented  in  Sections  2 and  3,  this  section  recom- 
mends some  constitutive  models  and  associated  failure  criteria  available 
for  practical  application  in  reinforced  concrete  structural  response 
prediction  by  the  finite  element  method.  Basically,  two  constitutive  models  are 
discussed  in  detail:  One,  proposed  by  Bazant  and  Bhat  (ref.  75)  and  the  other 
proposed  by  Chen  and  Chen  (refs.  70  and  71). 

1 . DISCRETE  MODELING 

Plain  concrete,  reinforcing  steel,  and  bond  link  relationships  will 
enable  detailed  representation  of  the  component  parts  of  reinforced  concrete. 

2.  PLAIN  CONCRETE 

The  Endochronic  Theory  of  Inelasticity  and  Failure  of  Concrete  proposed 
by  Bazant  and  Bhat  (ref.  75)  provides  the  most  complete  constitutive  law 
capable  of  response  prediction  for  varied  multistress  states.  Although 
the  theory  suitably  fits  finite  element  analysis,  incorporation  would 
require  considerable  modification  to  current  programs  available. 

Chen  and  Chen's  model,  though  less  accurate,  would  fit  well  in  standard 
finite  element  programs. 

a.  Constitutive  Relationship  Proposed  by  Bazant  and  Bhat: 

The  basic  stress-strain  relation  for  concrete  suggested  by  Bazant 
and  Bhat  is 


75 


1.  Deviatoric  Increment  component 


deij  = 


dS, . 

-^+deij" 


de  " = JJ  d 7 

e1j  2G  az 


05) 


2.  Volumetric  increment  component: 


= $ *de" 


de" 


dX  + ^ • — + d£° 
3K  t. 


(16) 


where 


= deviatoric  strain  and  stress  increments 

= volumetric  strain  and  stress  increments 

= deviatoric  stress  components 

= intrinsic  time.  Details  of  this  parameter  are 
in  paragraph  6. 

= inelastic  dilatancy;  inelastic  volumetric  strain. 

= stress-independent  inelastic  strain  (thermal 
dilatation  or  shrinkage) 

= updated  shear  and  bulk  moduli. 

= relaxation  time.  The  term  of  (dt/ij)  relates  to 
the  creep  effect. 

The  procedures  to  apply  the  stress-strain  relation  in  the  finite 
e lamest  method  suggested  is  as  follows: 

1.  After  the  linearization  problem  is  solved,  the  displacement 
increments  and  the  corresponding  strain  increments  in  the 
it  time  step  become  known:  Ar , Ac  ^ . 


"v  dsu 

de,  do 

S. . 
1J 

dz 

dX 

de° 

G,  K 

T. 


2.  Evaluate  the  parameter  A£  from  devlatorlc  strain  Increment  Aej 


6C  = /J2(Ae)  = A/2  Aeij-Ae^j 


3.  Compute  the  value  of  An  from  the  parameter  A £ by 


An  = F(e,a)A£ 


in  which 


F(e,a)  = a„  + 


V2 

a2[J2(e.)l  F3(maxe ) 
1 - aiMa)  - a3I 3(a) 


a0  = 0.2 


a3  = 0.4/f^ 
a 2 = 2500 
as  = 5/fl3 
n = n + An 
Co  = 0.0001 


4.  Calculate  Ac  from  An  by 


Ac  = An/f(n) 


= l + Bjn  + B2n2 


5.  Relaxation  time  tj  for  creep: 


It  may  be  noted  that  the  strain  Increment  calculated  in  each 
time  step  shall  be  decomposed  into  devlatoric  strain  components 
Ae.^  and  volumetric  strain  components  Ac. 

11.  Total  stress  increment: 

Aa  . . = AS..  + Aa  (33) 

' J • J 

b.  Constitutive  Equations  Suggested  by  Chen  and  Chen 

The  concrete  is  considered  an  elastic  plastic  strain-hardening 
and  fracture  material.  The  theory  uses  a failure  surface,  initial 
discontinuous  surface,  and  subsequent  discontinuous  surface  or  loading 
surface  to  define  the  behavior  (Figure  42).  Stress  states  causing 
complete  rupture  and  inability  to  resist  further  loading  define  the 
failure  surface.  The  initial  yield  surface  can  only  be  reached  by 
elastic  action;  unloading  within  this  surface  causes  no  permanent 
deformation.  Straining  beyond  this  surface  causes  permanent  deforma- 
tion and  modifies  the  shape  of  the  surface  at  which  subsequent  straining 
causes  permanent  deformation.  This  new  surface  is  the  loading  surface. 

Both  the  failure  and  initial  discontinuous  surface  depend  on  the 
devlatoric  stresses  and  hydrostatic  pressure  and  are  expressed  in  terms 
of  I i , and  J2: 


, . K2  K2  , 

' 3 J*  - 36  V * 


^ Ji  ‘ 36 


1 a T = r2 

3 U i T U 

(34) 

1 Vi  = x2o 

(35) 
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where  the  sign  of  the  third  term  Is  positive  In  compression  regions 
and  negative  In  tension-compression  regions,  and  AQ,  tq,  Ay,  and 
are  material  constants  which  can  be  determined  from  simple  tests. 

The  value  of  < is  set  to  3 based  on  experimental  results.  The  sub- 
sequent loading  surfaces  are  bounded  by  the  initial  discontinuous 
and  failure  surfaces,  and  assume  a similar  shape;  It  translates 
along  the  axis  and  expands  Isotropically: 


Using  the  normality  condition,  the  Incremental  stress-strain 
relations  are  derived  directly  from  the  loading  surface.  In  the  three- 
dimensional  case,  the  incremental  relation  is 


■ 


where 

l = { ( 1 - 2v ) (2n2  J2+3p2)  + 9vp2}  = _ _ ^l+v.)  ^~2v^  /(i^  j + 3ps) 

E 2 

$11  = {(1 -2v)  (nSxx  + p)  + 3vp}2 

$12  = {(l-2v)  (nSxx  + p)  + 3vp}  {(l-2v)  (nSyy  + p)  + 3vp} 

$13  = { ( 1 -2v)  (n$xx  + p)  + 3vp)  {(l-2v)  (nSzz  + p)  + 3vp} 

$14  = { ( 1 -2v)  (nSxx  + p)  + 3vp}  { (1  -2v)  riTxy} 

$15  = {(l-2v)  (nSxx  + p)  + 3vp}  { ( 1 -2v)  rrryz} 

$i6  = { ( 1 -2v)  (nSxx  + p)  + 3vp}  { (1-2v)  nxzx} 

$22  = { ( 1 -2v)  (nSyy  + p)  + 3vp}2 

$23  = {(l-2v)  (nSyy  + p)  + 3vp}  {(l-2v)  (nSzz  + p)  + 3vp}  (37) 

$24  = ( (l-2v)  (nsyy  + p)  + 3vp)  {(1-2v)  ntxy} 

$25  = {(l-2v)  (nSyy  + p)  + 3vp)  {(l-2v)  nxyz> 

$26  = { ( 1 - 2v ) (nSyy  + p)  + 3vp)  {(l-2v)  nxzx} 

$33  = ((l-2v)  (nS2Z  + p)  + 3vp}2 

$34  = {(l-2v)  (nSZ2  + p)  + 3vp)  { ( 1 - 2v ) nxxy} 

$35  = { ( 1 - 2v ) (nSzz  + p)  + 3vp}  { ( 1 - 2v ) nxyz} 

$36  = ((l-2v)  (nSZ2  + p)  + 3vp}  { ( 1 -2v ) nx2x} 

$44  = {(1-2v)  HXxy}2 
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$4  5 

- (O' 

-2v) 

nV 

{(1 

-2v) 

nxyz} 

$46 

= {(1 

-2v) 

nV 

(0 

-2v) 

$55 

= {(1- 

-2v) 

nTyz,! 

$56 

= {(1- 

■2v) 

nTyz} 

{(1 

-2v) 

^zx’ 

$66 

= {(1- 

-2v) 

nx  }2 
zx 

All  parameters  can  be  easily  calculated,  excepting  the  value  of  H. 
The  value  of  H,  defined  as  the  strain-hardening  rate,  can  be  written 
in  terms  of  the  slope  of  the  equivalent  stress-strain  curve: 

H = 2 a H'  (38) 

e 


where 

H1  = dae/dep 

ag  = the  equivalent  stress  = /f 

The  ofi  ~ ep  curves  in  the  compression  region  and  the  tension-compression 
region  are  different.  These  ae  ~ ep  curves  are  assumed  to  be  identical 
in  form  and  shall  be  obtained  from  tests.  The  curves  are  nonlinear  in  nature. 
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Since  the  derivation  of  the  incremental  stress-strain  equation  was 


based  on  the  elastic-plastic  model,  any  finite  element  program  for 
elastic-plastic  analysis  may  incorporate  this  material  model  without 
much  difficulty. 

Both  the  endochronic  theory  derived  by  Bazant  and  the  elastic- 
plastic,  strain  hardening  and  fracture  material  model  proposed  by  Chen 
and  Chen  have  been  presented.  Both  of  them  appear  to  be  promising  for 
finite  element  analysis.  However,  the  Bazant  model  provides  a better 
degree  of  fit  with  experimental  data  than  the  Chen  and  Chen  model.  The 
endochronic  theory  has  been  compared  extensively  to  different  experi- 
mental tests,  i.e., 

1.  Uniaxial  compression  test  by  Hognestad  et  al . (Fiqure  43). 

2.  Biaxial  stress-strain  test  by  Kupfer  et  al.  (Figures  44,  45,  46), 
and  biaxial  failure  envelope  (Figure  47). 

3.  Triaxial  tests  for  low  and  medium  strength  concrete  by  Balmer, 
1949,  and  Ri chart,  et  al.,  1928  (Figures  48  and  49). 

4.  Torsion-compression  tests  by  Bresler,  Pister  and  Gorde, 

Helony  (Figure  50). 

5.  Uniaxial  cyclic  compression  tests  by  Sinha  et  al.,  1964,  and 
Shak  et  al.,  1970  (Figure  51). 

6.  Creep  tests  (Figure  52,  Figure  53). 

Chen  and  Chen's  model  has  been  compared  only  to  the  biaxial  tests 
conducted  by  Kupfer  et  al.  (ref.  35).  The  figures  provided  by  Bazant 
and  Chen  and  Chen  are  given  herein  (Figures  42,  54,  55,  56).  Obviously, 
Bazant' s model  is  better  than  Chen  and  Chen's.  However,  it  shall  be  noted 
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that  the  derivation  of  Chen  and  Chen's  model  was  mainly  based  on  the 
classical  elastic-plastic  theory.  It  is  in  that  field  that  the  finite 
element  method  has  been  applied  most  successfully.  Further,  the 
incremental  stress-strain  relations  for  both  two-dimensional  and 
three-dimensional  stress  states  have  been  developed  in  a matrix  form 
that  is  easily  incorporated  into  any  existing  programs. 

3.  MATERIAL  MODEL  FOR  STEEL 


Compared  to  plain  concrete,  the  material  property  of  steel  is  very 
well-defined.  Current  yield  criterion  and  stress-strain  relations  for 
steel  have  been  demonstrated  to  be  adequate  in  predicting  steel  behavior 
even  in  the  plastic  range.  An  elastic-plastic  model  with  von  Mises  yield 
criterion  and  the  associated  flow  rule  may  be  assumed  for  reinforcing  steel. 
The  von  Mises  criterion  can  be  used  assuming  either  perfectly  plastic 
conditions  or  strain-hardening  situations.  The  procedures  and  formulations 
for  the  finite-element  method  applied  in  elastic-plastic  analysis  have  been 
well -developed.  It  would  not  be  difficult  to  handle  steel  in  reinforced 
concrete  analysis. 


The  reinforcing  steel,  in  finite  element  analysis,  may  be  treated  as 
individual  truss  bars,  or  assumed  to  be  distributed  uniformly  (smeared) 
over  the  entire  element  region  and  represented  by  2-D  or  3-D  solid  finite 
elements. 
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Strain- hardening  model 


Figure  42.  Chen  and  Chen's  Model 
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Volumetric  Strain,  ev  in  10-4 
Figure  45.  Volume  Change  in  Compression  Tests 
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Figure  47.  Biaxial  Failure  Envelope 
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Axial  Strain 


Figure  49.  Triaxial  Test  for  Medium  Strength  Concrete 


Torsion-Compression  Tests 


Cyclic  and  Static  Creep 
Whaley,  Neville,  1973 

a(t)  = aQ+  sin  wt 
aQ=0.35  f'c 


A a = f'  = 0.4 
c 


Endochronlc  theory 

f ■ 5650  psl  (39.0  MN/m2) 


Time  In  Hours 


Uniaxial  cyclic  compression 
low  stress 

(cyclic  creep)  /A 


V7~ 

/ A 


/s'/'  &/'/// 

•&  jP/'/ ,Z/v  vK' 


s/ySW'Y'' 

Hehmel,  Kern,  1962 

/' /sf  Endochronlc  theory 

.V  /Y''  f'  = 4860  psl  (33.51  MN/m2) 


Strain,  in  10 

Figure  52.  Cyclic  Creep  Tests 


Figure  53.  Nonlinear  Constant  Stress  Creep  and  Long  Time  Failure 


Figure  54.  Biaxial  Failure  Envelope 
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Library  Branch  G-043 
Washington,  DC  20545 

1 ATTN:  Doc  Con  for  Class  Tech  Lib 

US  Ehergy  Rsch  & Dev  Admin 
Nevada  Operations  Office 
P.0.  Box  14100,  Las  Vegas,  NV  89114 
1 ATTN:  Doc  Con  for  Tech  Lib 

Dept  of  the  Interior 
Bureau  of  Mines 

Bldg  20,  Denver  Federal  Center,  Denver,  CO  80225 
1 ATTN:  Tech  Lib 

Department  of  Defense  Contractors 

Aerospace  Corporation 
P.0.  Box  92957,  Los  Angeles,  CA  90009 
1 ATTN:  Tech  Info  Services 

Agbabian  Associates 

250  N.  Nash  St. , El  Segundo,  CA  90245 
1 ATTN:  M.  Agbabian 

Analytic  Services,  Inc. 

5613  Leesburg  Pike,  Falls  Church,  VA  22041 
1 ATTN:  George  Hesselbacher 
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Applied  Theory,  Inc. 

1010  Westwood  Blvd,  Los  Angeles,  CA  90024 
1 ATTN:  John  G.  Trulio 

Artec  Associates,  Inc. 

26046  Eden  Landing  Road , Hayward  CA  94545 
1 ATTN:  Steven  Gill 

Avco  Research  & Systems  Group 

201  Lowell  Street,  Wilmington  MA  01887 
1 ATTN:  Research  Lib  A830  Rm  7201 

Bat  telle  Memorial  Institute 

505  King  Ave. , Columbus  OH  43201 
1 ATTN:  Tech  Lib 

BDM  Corporation 

1920  Aline  Ave,.  Vienna  VA  22180 
1 ATTN:  Tech  Lib 

Ted  Belytschko 

6304  No.  Hiawatha  Ave. , Vienna  VA  22180 
1 ATTN:  Ted  Belytschko 

Boeing  Company 

P.0.  Box  3607,  Seattle,  WA  98124 
1 ATTN:  Aerospace  Lib 

Brown  Engineering  Co. , Inc 

Cummings  Research  Park,  Huntsville  AL  35807 
1 ATTN:  Manu  Patel 

California  Institute  of  Technology 

1201  E.  California  Blvd,  Pasadena  CA  91109 
1 ATTN:  Thomas  J.  Ahrens 

California  Research  & Technology,  Inc. 

6269  Variel  Ave. , Woodland  Hills  CA  91364 
1 ATTN:  Tech  Lib 

Cal span  Corporation 

P.0.  Box  235,  Buffalo  NY  14221 
1 ATTN:  Tech  Lib 
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Civil/Nuclear  Systems  Corp. 

1200  University  Blvd. , NE. , Albuquerque  NM  87102 
1 ATTN:  Tech  Lib 

Dayton,  University  of 
Industrial  Security  Super  KL-505 
300  College  Park  Ave. , Dayton,  OH  45409 
1 ATTN:  Hallock  F.  Swift 

University  of  Denver 
Colorado  Seminary 
Denver  Research  Institute 
p.O.  Box  10127,  Denver,  CO  80210 
1 ATTN:  Sec  Officer  for  Tech  Lib 

EG&G , Inc. 

Albuquerque  Division 
P.O.  Box  10218,  Albuquerque,  NM  87114 
1 ATTN:  Tech  Lib 

Electric  Power  Research  Institute 
3412  Hillview  Ave. , Palo  Alto,  CA  94303 
1 ATTN:  George  Sliter 

Engrg  Decision  Analysis  Co. , Inc. 

2400  Michelson  Dr. , Irvine,  CA  92715 
1 ATTN:  R.  P.  Kennedy 

Franklin  Institute 

20th  Street  and  Parkway,  Philadelphia,  PA  19103 
1 ATTN:  Zenons  Zudans 

General  Electric  Co.  , TEMPO-Center  for  Advanced  Studies 
816  State  St.,  (P.O.  Drawer  QQ),  Santa  Barbara,  CA  93102 
1 ATTN:  DASIAC 

General  Research  Corporation 
P.O.  Box  3587,  Santa  Barbara,  CA  93105 
1 ATTN:  Benjamin  Alexander 

H-Tech  Laboratories,  Lad. 

P.O.  Box  1686,  Santa  Monica,  CA  90406 
1 ATTN:  B.  Hartenbaum 
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I IT  Research  Institute 

10  West  35th  St.,  Chicago,  IL  60616 
1 ATTN:  Tech  Lib 

Institute  for  Defense  Analyses 

400  Army-Navy  Dr.,  Arlington,  VA  22202 
1 ATTN:  IDA  Librarian  Ruth  S.  Smith 

J.  H.  Wiggins,  co. , Inc. 

1650  S.  Pacific  Coast  Highway,  Redondo  Beach,  CA  90277 
1 ATTN:  John  Collins 

Kaman  Avidyne,  Division  of  Kaman  Sciences  Corp. 

83  Second  Ave. , Northwest  Industrial  park,  Burlington,  MA  01803 
1 ATTN:  Tech  Lib 

Kaman  Sciences  Corporation 

P.0.  Box  7463,  Colorado  Springs,  CO  80933 
1 ATTN:  Library 

Karagozian  and  Case 

6330  N.  Figueroa  St.,  Los  Angeles,  CA  90042 
1 ATTN:  John  Karagozian 

Lockheed  Missiles  & Space  Co. , Inc. 

P.0.  Box  504,  Sunnyvale,  CA  94088 
1 ATTN:  Tech  Lib 

Martin  Marietta  Aerospace,  Orlando  Div 

P.0.  Box  5836,  Orlando,  FL  32805 
1 ATTN:  G.  Fotieo 

McDonnell  Douglas  Corp. 

5301  Bolsa  Ave. , Huntington  Beach,  CA  92647 
1 ATTN:  Robert  W.  Halprin 

Merritt  Cases,  Inc. 

P.0.  Box  1206,  Redlands,  CA  92373 
1 ATTN:  Tech  Lib 

Newmark,  Nathan  M. , Consulting  Engineering  Services 

1211  Civil  Engineering  Bldg.,  Rm  B106A,  Univ  of  Illinois 


1 ATTN:  Nathan  M.  Newmark 
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TRW  Systems  Group 

One  Space  Park,  Redondo  Beach,  CA  90268 
1 ATTN:  Tech  Info  Center/S-1930 

1 ATTN:  Peter  K.  Dai,  Rl/2170 

1 ATTN:  Norm  Lipner 

TRW  Systems  Group,  San  Bernardino  Operations 
P.0.  Box  1310,  San  Bernardino,  CA  92402 
1 ATTN:  E.  Y.  Wong,  527/712 

TRW  Systems  Group 

Room  712,  Bldg  572, Norton  AFB,  CA  92409 
1 ATTN:  Gregory  D.  Eilcher 

Universal  Analytics,  Inc. 

7740  W.  Manchester  Bid. , Playa  Del  Rey,  CA  90291 
1 ATTN:  E.  I.  Field 

URS  Research  Co. 

155  Bovet  Rd. , San  Mateo,  CA  94402 

1 ATTN:  Tech  Lib 

Eric  H.  Wang  Civil  Engineering  Rsch  Fa,  University  Station 
Box  188,  University  of  NM,  Albuquerque,  NM  87121 

2 ATTN:  Jerry  Berglund 

Weidlinger  Assoc.  Consulting  Engineers 
110  East  59th  St. , New  York,  NY  10022 
1 ATTN:  Melvin  L.  Baron 

Weidlinger  Associates  Consulting  Engineers 
Suite  245,  3000  Sand  Hill  Road,  Menlo  Park,  CA  94025 
1 ATTN:  F.  S.  Wong 

1 ATTN:  J.  Isenberg 

Dept  of  Civil  Engineering 
University  of  California 
1 ATTN:  Prof  Alex.  C.  Scordelis 

729  Davis  Hall 
Berkeley,  CA  94720 

Purdue  University 
School  of  Civil  Engineering 
West  Lafayette,  IN  47907 
1 ATTN:  Prof  W.  F.  Chen 
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Pacifica  Technology 

P.0.  Box  148,  Del  Mar,  CA  92014 

ATTN:  G.  Kent 

ATTN:  Prof.  Robert  Dunham 

Physics  International  Co. 

2700  Merced  St. , San  Leandro,  CA  94577 
ATTN:  Doc  Con  for  Tech  Lib 


R & D Associates 

P.0.  Box  9695,  Marine  Del  Rey,  CA  90291 
ATTN:  Tech  Lib 


Rand  Corporation 

1700  Main  St. , Santa  Monica,  CA  90406 
ATTN:  Tech  Lib 

Science  Applications,  Inc. 

2201  San  Pedro  NE,  Albuquerque,  NM  87110 
ATTN:  J.  L.  Bratton 

Science  Applications,  Inc. 

P.  0.  Box  2351,  La  Jolla,  CA  92038 
ATTN:  Tech  Lib 

Southwest  Research  Institute 

P.0.  Box  Drawer  28510,  San  Antonio,  TX  78284 

ATTN:  Wilfred  E.  Baker 


Stanford  Research  Institute 

333  Ravenswood  Ave. , Menlo  Park,  CA  94026 

ATTN:  Carl  Peterson 

Systems,  Science  and  Software,  Inc. 
p.O.  Box  1620,  La  Jolla,  CA  92038 
ATTN:  Tech  Lib 
ATTN:  G.  Hegemeir,  Room  2091 

Terra  Tek,  Inc. 

420  Wakara  Way,  Salt  Lake  City,  UT  87108 
ATTN:  Tech  Lib 

Tetra  Tech,  Inc 

630  North  Rosemead  Blvd. , Pasadena,  CA  911u/ 
ATTN:  Tech  Lib 
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PMB  Systems  Big 
500  Sansome  Street 
Ran  Prancisco  94111 


